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We performed fluctuation analysis by means of the local scaling dimension for the strength function 
of the isoscalar (IS) giant quadrupole resonance (GQR) in ^"^Pb where the strength function is 
obtained by the shell model calculation including Iplh and 2p2h configurations. It is found that at 
almost all energy scales, fluctuation of the strength function obeys the Gaussian orthogonal ensemble 
(GOE) random matrix theory limit. This is contrasted with the results for the GQR in '^''Ca, where 
at the intermediate energy scale about 1.7 MeV a deviation from the GOE limit was detected. It is 
found that the physical origin for this different behavior of the local scaling dimension is ascribed 
to the difference in the properties of the damping process. 
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I. INTRODUCTION 

Giant resonances, excited by various probes, show, at 
an initial stage of the excitation process, a regular motion 
with a definite vibrational frequency [Ij, |2i] . These regular 
motions are then damped due to the coupling with a huge 
number of background states, and finally the so called 
compound states are realized. 

We now have understood the both ends of these pro- 
cesses: The frequency of the giant resonance, for in- 
stance, can be well evaluated by the random phase 
approximation (RPA). Compound states, on the other 
hand, are also well described by the random matrix the- 
ory with the Gaussian orthogonal ensemble (GOE) 
which characterizes a classical chaotic motion. 

It is still not well understood, however, how the dy- 
namics changes from regular to chaotic '5j. In order to 
answer this question, it is very useful to study the fluctu- 
ation properties of the strength functions: The structure 
at the large energy scale of the strength function cor- 
responds to the behavior of the initial stage, while the 
fluctuation properties at small energy scale correspond 
to the long time behavior. 

We proposed and have used a novel fluctuation analy- 
sis based on the quantity we call the local scaling dimen- 
sion to study the fluctuation properties of the strength 
functions [j^]. This method is devised to quantitatively 
characterize the fluctuation at each energy scale, and is 
suitable for the investigation of the fine structure of the 
strength function. 

The strength distribution of giant resonances and its 
fluctuation have also been studied experimentally. Re- 



cently, the fine structure of the strength distribution of 
the giant quadrupole resonance (GQR) in ^°^Pb M, H, Q 
or the Gamow- Teller resonance (GTR) in ^°Zr were 
measured and theoretical analysis has also been done. 

In the previous paper we investigated the GQR 
in ""^Ca, where the strength function was calculated 
by means of the second Tamm-DancofF approximation 
(TDA), namely, the Iplh and 2p2h model space is in- 
cluded. The results of the local scaling dimension anal- 
ysis were as follows: At small energy scale, the behavior 
of the local scaling dimension is almost the same as that 
of the GOE, which exhibits the complexity of 2p2h back- 
ground states. On the other hand, a clear deviation from 
the GOE was found at the intermediate energy scale and 
it was found that this energy corresponds to the spread- 
ing width of Iplh states. Hence, we can say that the 
spreading width of Iplh states is detected as deviation 
from the GOE limit in ^OCa. 

For ^"Ca the Landau damping is important for the 
damping process of the giant resonance. Namely, the 
strength is first fragmented over a wide range of Iplh 
states, and this fragmentation characterizes a global pro- 
file of the total strength function. 

However, as the mass of nuclei increases, the relative 
importance of the Landau damping may change. Accord- 
ingly, 2p2h states may also contribute to the global profile 
of the strength function. Therefore, it is very important 
to investigate how the difference between the damping 
process of light nuclei and that of heavy nuclei does affect 
the properties of the ffuctuation of the strength function. 

In this paper, we study the isoscalar (IS) GQR of 
^"^Pb, where the strength function is calculated with the 
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second TDA in the same manner as in ''^'^Ca, and study 
the fluctuation of the strength function by means of the 
local scaling dimension. Comparing results with those 
of ^"Ca we would like to clarify which properties of the 
damping process are reflected in the fluctuation of the 
strength function and make clear the physical origin of 
the difference. 

This paper is organized as follows: In Sec. [Ill we briefly 
explain the local scaling dimension. The strength func- 
tion for IS GQR in 2°*Pb is calculated in Sec. nil where 
the adopted Hamiltonian and the model space are shown. 
In Sec. IIVI we discuss the nearest-neighbor level spacing 
distribution, A3 statistics as well as a histogram of the 
strength distribution. In Sec. |Vl we apply the local scal- 
ing dimension to the IS GQR strength function in ^°*Pb. 
Detail of damping process is studied in Sec. lVIi where the 
physical origin for the difference of the fluctuation prop- 
erty of the strength function between '*°Ca and ^"^Pb is 
also discussed. Finally, Sec. lVIIl is devoted to conclusion. 
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FIG. 1: Calculated strength function of the IS GQR in^°**Pb. 
Dotted curve shows the smooth strength function by means 
of the Strutinsky method with the smoothing width 0.2 MeV. 



II. LOCAL SCALING DIMENSION 

We briefly explain the local scaling dimension. 
Refs. for details. 

The strength function is expressed as [l2l | 



See 



(1) 



Here Ei and Si denote the energy and the strength of 
exciting the ith energy level, respectively. Strengths are 
normalized as ^ ■ 5^ — 1. 

To study the fluctuation at each energy scale, we con- 
sider binned distribution of the strength by dividing 
whole energy interval under consideration into N bins 
with length e. Strength contained in nth bin is denoted 
by Pn, 
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Pn = > Si 
iGnth bin 



(2) 



To characterize the distribution of the binned strengths, 
we introduce the moments of p„ , which are called in lit- 
erature the partition function Xm(e) defined by 



Xm(e) 



N 



71 = 1 



(3) 



Finally, by extending the idea of the generalized fractal 
dimensions [13, [3] to non-scaling cases in a straightfor- 
ward way, we can define the local scaling dimension as, 



— 1 d log e 



(4) 



Since the local scaling dimension has a definite physical 
meaning similar to that of the generalized fractal dimen- 
sion, the value of can quantitatively characterize 



the fluctuation of the strength function at each energy 
scale e. 

In the actual calculation of the local scaling dimension, 
we define it by means of the finite difference under the 
change of a factor 2, 

n / /o N 1 logXm(2e) - logxm(e) 

Dm(V2e) = — , 5) 

m — 1 log 2 

rather than the derivative in Eq. 



III. NUMERICAL CALCULATION OF 
STRENGTH FUNCTION 

We calculated the strength function of the IS GQR in 
208pb within the second TDA including the Iplh and 
2p2h excitations. Single-particle wave-functions and en- 
ergies were obtained for a Woods-Saxon potential includ- 
ing the Coulomb interaction. The effective mass parame- 
ter m*/m, which scales the Woods-Saxon single-particle 
energies ews as Ehf = £ws / {m* / ra) to simulate the bare 
(Hartree-Fock) single-particle energies Ehf, is set to be 1 
in this calculation. 

As the residual interaction, the Landau-Migdal-type 
interaction including the density-dependence was 
adopted. The model space was constructed in terms of 
single-particle states within the four major shells, two 
below and two above the Fermi surface, and included all 
Iplh states and 2p2h states whose unperturbed energies 
are less than 15MeV. Resultant number of Iplh states 
and 2p2h states are 39 and 8032, respectively. We diago- 
nalized the Hamiltonian within this model space and ob- 
tained the strength function for the isoscalar quadrupole 
operator. 

Figure [T] shows the calculated strength function. The 
average of the excitation energy weighted by the strength 
is about 10.5 MeV, and the standard deviation around 
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FIG. 2: The nearest-neighbor level spacing distribution for 
(a) ^°Ca and (b) ^^^Pb. For ^"^Ph 3321 levels between 9.9 
MeV and 13.1 MeV, while for *°Ca 804 levels between 20 MeV 
and 30 MeV are considered. Level spacings were unfolded by 
the Strutinsky method with a smoothing width 0.5 MeV for 
^°®Pb and 5.0 MeV for ''"Ca, respectively. The solid curve 
represents the Wigner distribution. 



FIG. 3: The histogram of the square-root of normalized 
strengths Sj^^ associated with IS GQR in (a) *"Ca and (b) 
^°*Pb. The solid curve represents the Porter- Thomas distri- 
bution which becomes a Gaussian when plotted as a function 
of Sj^'^ . See the caption of Fig. [2] for the number of consid- 
ered levels and also see Sec. IV Al for the normalization of the 
strengths. 



the average is about 2.6 MeV, where all levels are consid- 
ered. The peak position lies at the same value as the av- 
erage. These values are consistent with the (p,p') exper- 
imental data Moreover, the agreement of the global 
shape with the experimental data is also good. The dot- 
ted curve in Fig. [1] represents the smooth strength func- 
tion by means of the Strutinsky method with the 
smoothing width 0.2 MeV. The value of the FWHM of 
this smooth strength function is 0.63 MeV. In order to 
quantitatively characterize the spreading of the strength 
function around the largest peak, the FWHM is more 
appropriate than the standard deviation |17| . Thus, we 
use the FWHM as a measure of the total width F of the 
strength function, which gives F = 0.63 MeV. 

Hereafter, when we estimate the value of the FWHM, 
the same procedure as above is adopted, namely, we 
calculate the FWHM for the smooth strength function 
by means of the Strutinsky method with the smoothing 
width 0.2 MeV. 
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FIG. 4: The A3 statistics for (a) ""Ca and (b) ^"^^Pb. The 
horizontal axis L shows the value of the energy interval for 
the unfolded spectrum. The solid curve represents the A3 for 
the GOE level fluctuation. See Fig. |5] for other parameters. 



IV. FLUCTUATION AT SMALL SCALE 



Before going to the detailed discussion of the local 
scaling dimension, we briefly show the results for other 
fluctuation measures: the nearest-neighbor level spacing 
distribution (NND), the strength distribution, and A3 
statistics. Here, the NND and the strength distribution 
are measures characterizing the fluctuation at small en- 
ergy scale limit. We present the results of ^°Ca as well 
as those of ^''^Pb for the sake of comparison. 

Figure [2] shows the NND. For both nuclei the NND 
follows the Wigner distribution well. We present the 
strength distribution in Fig. [3] where a histogram of the 
square-root of normalized strengths is plotted. We also 



flnd that for both "''Pb and °Ca the distribution follows 
the Porter-Thomas one rather well. These two figures in- 
dicate that for both nuclei the fluctuation of the strength 
as well as that of the energy level spacing is governed by 
the GOE at least at small energy scale limit as expected. 

Figure [5] shows the A3 statistics. We again find that at 
small energy range the A3 follows the GOE line for both 
^"^Pb and ^"^Ca, although at intermediate energy scales, 
Lmax ^ 20 or 15 for ^ospb or "^Ca, respectively, the A3 
starts to deviate from the GOE line to upward. 
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RESULTS OF LOCAL SCALING DIMENSION 



A. Calculation of the local scaling dimension 

Since we are not interested in the global shape of 
the strength function, we actually adopt the normalized 
strength function S(E) for the fluctuation analysis as in 
the case of '^''Ca pjj. The normalized strength function 
S{E) is given by 



(6) 



where the normalized strength Si of the ith level is de- 
fined by 



(7) 



Here, p{E) and S{E) denote the level density and the 
strength function, respectively, smoothed by the Strutin- 
sky method ^16,] . A/" is a normalization factor to guarantee 

We determine the width parameter uj of the Struti- 
nsky smoothing function as follows: We note that the 
smoothed strength function S{E) should represent the 
global profile of the original strength function S{E) at 
large energy scale, but at the same time, we would like 
to choose UJ as large as possible since we do not want to 
wash out the fluctuations at smaller energy scales. Fig- 
ure [H] shows the FWHM of the smoothed strength func- 
tion S{E) / p{E) as the function of the smoothing width 
Lo. The linear increase of the FWHM at large values of 
UJ > 0.6 MeV indicates that the value of uj is too large, 
while with smaller values uj < 0.5 MeV the FWHM stays 
at an approximately constant value, reflecting the total 
width. We therefore adopt 0.5 MeV as the value of the 
smoothing width uj in order to satisfy the above require- 
ments. 

We use the equidistant energy level Ei in Eq. 
namely, Ei = id, where d denotes the average level spac- 
ing. Finally, we adopted the energy range from 9.9 MeV 
to 13.1 MeV, where 3321 levels are included. 

The normalized strength function is plotted in Fig. [B] 
The local scaling dimension is derived from this normal- 
ized strength function. 



B. Behavior of the local scaling dimension 

Figure [7] (a) and (b) represent the partition function 
and the local scaling dimension, respectively, of IS GQR 
in ^''^Pb. The horizontal axes in both figures represent 
the bin width e of energy in unit of d, where d rep- 
resents the average level spacing over the energy range 
9.9 - 13.1 MeV {d = 0.96 keV). The partition function 
clearly deviates from the linear relation in the log-log 
plot. This means that for the GQR strength function 
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FIG. 5: FWHM of the smoothed strength function 
S{E)/p{E) of IS GQR in 2°*Pb as a function of smoothing 
width UJ used in the Strutinsky method. The dotted line is 
fitted to data and gives ~ 0.98a; + 0.2. 
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FIG. 6: Normahzed strength function Eq. (|6j of IS GQR 
in^°*Pb. Smoothing width uj = 0.5 MeV was used. 



the self-similar property does not hold. We can also see 
a more detailed structure in the figure of the local scaling 
dimension. At the smallest energy scale e ~ d, the value 
of the local scaling dimension is small, D2 — 0.35, which 
means that the fluctuation is very large at small energy 
scales. As the energy scale or the bin width increases, 
the values of D„i{e) monotonically increase. Finally, at 
about e ~ lOOd the values of £'m(e) converges to unity, 
which indicates that at large energy scales, the strength 
function appears smooth. The most important feature in 
Fig.[7](b) is that the local scaling dimension for ^°*Pb al- 
most follows the GOE line at almost all the energy scales. 

This should be contrasted with the case of '*°Ca [Uj : 
The partition function and the local scaling dimension 
for '"'Ca are shown in Fig. [7] (c) and (d), respectively, 
for a comparison. When the energy scale is small, the 
local scaling dimension almost follows the GOE line. As 
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FIG. 7: Partition function (a) and local scaling dimension (b) 
for the IS GQR in '^°^Ph, and those in ^°Ca are also shown 
at (c) and (d). Curves in each figure correspond to m = 2 - 5 
from upper to lower. Dotted curves in (b) and (d) represent 
D2(e) for the GOE. 

the energy scale increases, however, we can find a dip 
and a deviation from the GOE line at about 1.7 MeV 
(Note that d = 12 keV for ''°Ca). We verified that an 
occurrence of the dip is not due to a statistical error. 
Moreover, further studies indicate that the energy where 
the minimum is located is approximately related to the 
value of the spreading width of Iplh states. 

Note that if we look only at the small energy scale 
limit or large energy scale limit, we can not find the dif- 
ference between ^'^^Pb and '"^Ca. Studies of fluctuation 
at intermediate energy scales lead to the finding of the 
difference. In the following we shall investigate the mech- 
anism which brings about the difference in fluctuations 
at intermediate energy scales. 

VI. STUDIES OF DAMPING PROCESS 

Let us now investigate origins of the difference between 
the cases of ""^Ca and ^°^Pb. In our previous study of 
the GQR in ^"Ca, we have shown that the behavior of 
the local scaling dimension, shown in Fig. [7] (d), can be 
interpreted in terms of the doorway damping mechanism. 
We here employ the same picture in order to clarify the 
damping mechanism of the GQR in ^°*Pb. 

The doorway damping mechanism consists of a two- 
step process which is illustrated in Fig. [5] The giant res- 
onance is spread over the Iplh states due to the Landau 




FIG. 8: Schematic drawing of the doorway damping mecha- 
nism of the giant resonance, and related quantities. 



damping, the width of which is denoted by Fl. The aver- 
age spacing of Iplh states is denoted by fipih. The Iplh 
states are considered here as the "doorway" states of the 
damping process. The Iplh states then couple to more 
complicated background states (2p2h states) through the 
residual two-body interaction. The coupling causes the 
spreading width of Iplh states, which we denote 712. We 
define the GQR TD state as the Tamm-Dancoff (TD) 
state with the largest quadrupole strength among all TD 
states, where the TD states mean the states obtained in 
the TDA, i.e., by the diagonalization within the model 
space limited to the Iplh configurations. The GQR TD 
state also couples to 2p2h states, and hence it should 
have the spreading width due to the coupling. This is 
similar to 712, but we introduce a separate symbol r2 
since the GQR TD state is a special state consisting of a 
coherent superposition of many unperturbed Iplh exci- 
tations. d2p2h is the average spacing of background 2p2h 
states. The residual interaction also acts among the 2p2h 
states, and the mixing among the 2p2h states causes a 
spreading width of the 2p2h states, which we denote 722. 

In the following we shall evaluate all these quantities 
in order to clarify the damping mechanism of the GQR in 
208pb (Sec. EE] and Sec. Ell]). We also study whether 
there are specific states among 2p2h states which strongly 
couple with the GQR mode (Sec. IVICj) and then discuss 
the difference of the nature associated with the fluctua- 
tion of strength function between ^"Ca and ^°^Pb (Sec 

Emj. 



A. Mechanism producing the total width 

1. Landau damping 

For ^"Ca, the Landau damping is important, so that 
the strengths are already fragmented in the Iplh levels. 
Therefore we first would like to investigate in ^°^Pb, how 
the strength is distributed in the TDA where only the 
Iplh states are included. 

Figure|5]shows the TDA strength function, which is ob- 
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FIG. 9: TDA strength function for the IS quadrupole oper- 
ator in ^°*Pb. See Fig.[T]for the dotted curve. 



FIG. 10: Strength function by neglecting all TD states except 
the GQR TD state. 3342 2p2h states lying in 9 MeV - 13 MeV 
are considered. See Fig. [l] for the dotted curve. 



tained by means of the TDA, namely by neglecting 2p2h 
states, of the IS quadrupole operator. Different from the 
case of ''"Ca, strengths in the GQR region is consider- 
ably concentrated on the single peak located at about 
10.7 MeV. Because of this, the TDA strength function is 
very different from the full strength function in Fig. [TJ 
At the same time, we also see only a small effect of the 
Landau damping. In fact, the strength concentration on 
the single peak at E = 10.7 MeV is 59% of the strengths 
in the energy interval 9-13 MeV. The Landau damp- 
ing width Tl may be evaluated in terms of a smoothed 
profile of the strength function plotted with the dotted 
curve in FigM Its FWHM reads 0.21 MeV. On the other 
hand, if we closely look at Fig. [HI we find that there is the 
second largest peak just below the largest one and that 
these two levels dominate the whole structure. The level 
spacing between these two levels can be considered as a 
typical spreading of strength and may be a more direct 
quantitative measure of the Landau damping width Fl: 
The level spacing 0.18 MeV gives Fl = 0.18 MeV. 



2. damping due to 2p2h states 

The Landau damping width Fl = 0.18 MeV is not 
enough to explain the total width F = 0.63 MeV of Sec. 
mil Then, we would like to study a role of 2p2h states in 
the damping process, namely, the fragmentation of the 
GQR TD state located at E ^ 10.7 MeV in Fig. [9] over 
2p2h states. We shall investigate the damping width F2 
caused by the coupling to 2p2h states. To estimate this 
width, we perform a calculation where we include only 
the GQR TD state and 2p2h states, where the coupling 
between the GQR TD state and 2p2h states as well as the 
interaction among 2p2h states are taken into account. 

Figure \W\ shows the resulting strength function. The 
estimated FWHM is 0.41 MeV, i.e., F2 = 0.41 MeV. 

If the Landau damping and the 2p2h damping are inde- 



pendent of each other, and neighboring TD states around 
the GQR TD states also have the same spreading width 
as F2, the following approximate relation holds: 



F : 



(8) 



The values, Fl = 0.18 MeV and F2 = 0.41 MeV, esti- 
mated above indeed satisfy this relation. Consequently, 
the total width F = 0.63 MeV is approximately explained 
as a sum of the Landau damping width Fl and the 2p2h 
damping width F2. 

The importance of the 2p2h damping is contrasted 
with the case of ^"Ca, where the total width can be ex- 
plained essentially by the Landau damping width, i.e., 
L-Fl. 



B. Spreading width of Iplh states and 2p2h states 

For the case of *°Ca, the strength is fragmented over 
many Iplh states by the Landau damping, and strength 
in each Iplh state is further spread due to the coupling 
with 2p2h states. Let us evaluate the spreading width 
712 of the Iplh states due to this couphng. We shall also 
evaluate the spreading width 722 of 2p2h states, which is 
caused by the residual coupling among 2p2h states. 

We evaluate 712 by using the strength functions of TD 
states as in Ref. ll|. Namely, we calculate the strength 
function of each TD state. Averaging the strength func- 
tions over whole TD states, we obtain Fig. [Tl] (a). The 
FWHM of this averaged strength function gives an eval- 
uation of the spreading width 712. We read 712 — 0.38 
MeV. (Note that we define 712 as the spreading width of 
TD states instead of that of unperturbed Iplh states.) 
The value of spreading width of 2p2h states 722 is also 
evaluated in the same manner. From Fig. [TT] (b) we also 
obtain 722 = 0.75 MeV as the estimate of the spreading 
width of 2p2h states. These results will be used in Sec. 

mm 
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FIG. 11: Averaged strength function of (a) TD states and 
(b) 2p2h states. Average was performed over levels lying in 9 

MeV - 13 MeV. The number of levels is 12 and 3342 for TD ^^^jg^^: TDA strength function for the IS octupole operator 
states and 2p2h states, respectively. 



For the sake of comparison, let us estimate 
the spreading width by assuming the Fermi golden 
rule. The root mean square of matrix elements be- 
t ween Iplh state s and 2p2h states is calculated as 
((lplh|V^i2|2p2h)2)i/2 ^ 9,3 X 10-3 MeV. Similarly, we 

calculate ((2p2h|\/22|2p'2h') 2)^/2 = 1.0 x 10"^ MeV. 
Since the level spacing of 2p2h states is c?2p2h = 1-2 
keV, the spreading widths 712 and 722 are approxi- 
mately estimated in the Fermi golden rule as 7^2^ = 
27r(lplh|V"i2|2p2h)2/d2p2h = 0.46 MeV and 7^3° = 
27r(2p2h|X^22|2p'2h')Vc^2p2h = 0.53 MeV, respectively, 
which arc in approximate agreement with the direct eval- 
uation within 30%. 



C. Search for strongly coupled states in 2p2h states 

In the picture of Fig.[5]2p2h states are assumed to play 
a role as the chaotic background and provide the GOE 
fluctuation to the strength function. However, if the 
GQR TD state couples with not all 2p2h states equally 
but specific states in 2p2h states strongly, there is a pos- 
sibility for this hierarchical structure in 2p2h states to 
give rise to a deviation from the GOE fluctuation. We, 
here, would like to investigate whether whole 2p2h states 
are rather equally coupled with the GQR TD state or 
whether there are specific states in 2p2h states which 
strongly couple with that state. 

As a candidate of such specific states, we can consider 
the low-en ergy surface vibration plus Iplh states: In 
Refs. [H [il Hi, HO, the importance of the cou- 
pling to the surface vibration in the wide range of damp- 
ing phenomena including the damping of a single parti- 
cle motion as well as that of giant resonances was dis- 
cussed. As for the giant resonance, which is composed 
of a coherent superposition of Iplh states, this means 
that the damping occurs via the coupling with the spe- 
cific 2p2h states, namely, the surface vibration plus Iplh 



(s.v.+lplh) states. 

Since our model does not assume the particle- vibration 
coupling a priori, it is not trivial whether our model 
also has a mechanism that enhances the coupling with 
the low-energy surface vibration. Therefore, we would 
like to study whether the s.v.-|-lplh states are particu- 
larly strongly coupled with the GQR TD state within 
our model. To do so, we calculate the FWHM of the 
following approximate strength function: 

^ ' n y ^S-c^„ +1722/2 J 

(9) 

where, Ec and uja denote the energy of the GQR TD state 
and the energy of the ath s.v.-|-lplh state, respectively. 
Vca represents the coupling matrix element between the 
GQR TD state and the s.v.-f Iplh state a. 

Only J'^ =2+, 3~ modes are included as surface vibra- 
tions: We took only the lowest TD state as J'^ = 2+ sur- 
face vibrational mode. On the other hand, we must pay 
attention to the collectivity of the octupole mode. Figure 
[T^ shows the TDA strength function for the IS octupole 
operator. Compared with the experimental data [23 |. the 
energy of the lowest state is too high, and strengths are 
fragmented over several states. Thus, we took into ac- 
count the lowest nine states for the octupole mode. Note 
that s.v.-|-lplh states thus defined are not orthogonal. In 
this sense Eq. ^ is an approximation which neglects the 
non-orthogonality. 

The strength function based on Eq. Q is presented in 

Fig. [H The width r^"'"-^ estimated by the FWHM is 
0.074 MeV. This value is significantly smaller than the 
width r2 = 0.41 MeV of the GQR TD state caused by 
the coupling to the whole 2p2h states. 

From the estimate by the Fermi golden rule, we can 
give more detailed comparison between the width for the 
case of s.v.-flplh states and that for the whole 2p2h 
states. It is noted in Table U that the spreading width 
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FIG. 13: Strength function of the GQR TD state evaluated 
by considering only surface vibration plus Iplh (s.v.+lplh) 
states based on Eq. 722 = 0.75 MeV is used. 



TABLE I: Averaged value of squared coupling matrix el- 
ements between the GQR TD state and surface vibra- 
tion plus Iplh states or the whole 2p2h states(third column), 
the associated spreading width Ff*' of the GQR TD state 
evaluated by the Fermi golden rule (fourth column), and the 
spreading width F2 estimated by the FWHM of the strength 
function based on Eq. @ (fifth column). Second column 
shows the number of states considered. The second row shows 
the results obtained by including only the s.v.-l-lplh states 
while the third row shows those for the case of the whole 
2p2h states. 



# (MeV^) 


F^° (MeV) 


F2 (MeV) 


s.v.-Hlplh 909 0.65 x 10"'' 


0.092 


0.074 


2p2h 3342 0.72 x 10"" 


0.38 


0.41 


r^'-''-^ = 0.074 MeV and Fa 


= 0.41 MeV 


are well ac- 


counted for by the estimate. 


In the Fermi 


golden rule 



the spreading width is governed by two factors; 1) the 
average value of squared couphng matrix elements V^^ 
between the GQR TD state and the states that couple 
to it, and 2) the level density of the coupling states. 
From Table [H we see that the large difference between 
the two widths simply reflects the difference between 
the number of s. v. -(-Iplh states 909 and 2p2h states 
3142 whereas the coupling strength of s.v.-|-lplh states 
Kq = 0-65 X 10"^ MeV^ is comparable to the couphng 
strength = 0.72 x 10"^ MeV^ for the whole 2p2h 
states. 

Table H] and Fig. [T3] suggest that our model does not 
contain the enhancement of the coupling with the sur- 
face vibrations in the damping of the GQR. Therefore we 
consider in the following the 2p2h states as background 
states which do not have specific structures. 



TABLE II: Values of physical quantities related to the damp- 
ing of the GQR for ""Ca and ^"^Ph. Unit of the energy is keV 
for all cases. 



r 


Fl 


F2 


712 




722 


d2p2h 


■*°Ca 4000 


4000 


1500 


1500 


500 


5200 


11 


^°*Pb 630 


180 


410 


380 


230 


750 


1.2 



D. Physical origin of the difference between Ca 



In the above subsections, we have evaluated the physi- 
cal quantities such as the various spreading widths, with 
which we have discussed the damping process of ^°Ca and 
^'^^Pb, especially the mechanism of producing the total 
width of the strength function. Here, using these quan- 
tities we would like to discuss the physical origin of the 
difference between the fluctuation of the strength fluctu- 
ation of '"^Ca and that of ^"^Pb. Table [Til summarizes 
the values of the above physical quantities related to the 
initial stage of the damping process for both ^"Ga and 

208pb. 

We have shown in our previous study 6] that the 
damping process through the doorway states causes large 
fluctuations which have characteristic energy scales, and 
that the fluctuations emerge in the local scaling dimen- 
sion. For instance, the energy scale of the spreading 
width 712 of the doorway states is the quantity which 
shows up prior to the other quantities. It is noted, how- 
ever, the size of the fluctuations depends on the mutual 
relations among the quantities mentioned above, and in- 
deed we have examined in the relations which are 
needed to detect the effect of the spreading width 712. 

It is trivial that the local scaling dimension can detect 
the spreading width when the spreading of Iplh states 
does not cause the overlap of these states, namely when 
712 < i^ipih. In addition to this case, the local scaling 
dimension still keeps the information of the spreading 
width even if the Iplh states start to overlap with each 
other, i.e. 712 ~ i'lpih- Studying more quantitatively 
with the use of the doorway damping model of Ref . Q , we 
found the condition to detect the effect of the spreading 
width as 

(A) 712 < 4i:>ipih. 

Furthermore, we need the second condition: 

(B) 712 < Fl. 

This simply means that the spreading width 712 of the 
doorway states (Iplh states) need to be smaller than the 
total width F. Since F ~ Fl + F2 and 712 ~ F2, the 
requirement 712 < F can be written as (B). In addition 
to (A) and (B), we need the third condition: 

(G) i^ipih < Pl. 
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FIG. 14: Schematic picture of the initial stage of the damping 
process for GQR in ^"^'Pb. 

This is because we need more than one doorway states 
within the the energy interval Fl in order to have fluc- 
tuating behavior in the strength function. 

Let us first look at the case of ''"Ca. From Table [III 
the relation 712 = S.Ofipih is derived, and this relation 
fulfills the condition (A). On the other hand, relations 
712 — O.SSFl and -Dipih — O.ISFl are also derived from 
Table HIl and these relations satisfy both conditions (B) 
and (C). As a result, in the case of ''"Ca, we can see a 
deviation from the GOE fluctuation in the local scaling 
dimension and indeed the energy scale where the devia- 
tion is seen is related to the value of 712. 

For ^°^Pb, on the other hand, we find in Table [Til that 
712 = 1.7Dipih, while F^ is smaller than 712 and I?ipih, 
i-e., 712 = 2.1Fl and I^ipih — 1.3Fl. The first relation 
satisfies the condition (A). The latter two relations, how- 
ever, break the condition (B) and (C). Accordingly, for 
the case of ^"*Pb, the deviation from the GOE due to 
the effect of 712 can not be seen. The situation in ^"^Pb 
is illustrated in FigdU Essential physical origin of this 
difference is that for ^°*Pb the Landau damping width 
is small compared with that of '*°Ca. The smallness or 
largeness of the value of the Landau damping width af- 
fects the fluctuation property of the strength function. 

VII. CONCLUSION 

We studied the fluctuation properties of the strength 
function of IS GQR for ^°^Pb by means of the local scal- 



ing dimension, and compared the results with those of 
""^Ca. The strength function was obtained by the second 
TDA including 2p2h states as well as Iplh states. For 
''"Ca, we find a fluctuation different from GOE around 
the energy scale which is approximately related to the 
spreading width of the Iplh states. On the other hand, 
for ^°^Pb we can not find the fluctuation different from 
the GOE at almost all the energy scales. 

The different behavior of the fluctuation detected by 
the local scaling dimension analysis is due to the differ- 
ence of the ratio of the Landau damping width Fl to the 
spreading width of the Iplh states 712. 

Recently, the analysis of the strength function of the 
IS GQR in ^°*Pb obtained by (p,p') inelastic scatter- 
ing experiment was performed by means of the wavelet 
transform Q. The authors suggest from the positions 
of the local maxima in the wavelet power that there ex- 
ist three energy scales in the fluctuation of the strength 
function: I. 120 keV, II. 440, 850 keV, III. 1500 keV. 
Existence of higher two energy scales is not inconsistent 
with our results, since our analysis says nothing about 
the fluctuation at about energy scale II, which may cor- 
respond to the total width F in our model, or higher 
energy scales. However, the existence of the smallest en- 
ergy scale ~ 120 keV may conflict with our results: If 
there is such an energy scale in our strength function, 
our analysis must detect it as a deviation from the GOE 
fluctuation. Therefore, it is very important to study the 
origin of this discrepancy. In particular, it is interesting 
to clarify the relation between two method, namely, the 
local scaling dimension and the wavelet power. Studies 
in this direction are now in progress. 
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